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We study the single-particle density of states of one-dimensional and two-dimensional quantum 
disordered systems with long-range interactions. We consider a l/y/r interaction in one dimension 
and a Coulomb interaction in two dimensions, which produce linear gaps in the density of states in 
both cases. We focus on the strong localization regime where the localization length is small but non- 
zero. We use an exact diagonalization technique for small system sizes and a perturbation approach 
for larger sizes. We find that, with both methods, the inclusion of a finite hopping contribution does 
not change the linear character of the gap, but reduces its slope, widening the gap. 
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I. INTRODUCTION 

Insulators can exhibit interesting transport properties. 
At very low temperatures localized electrons jump via 
phonon-assisted hopping and the transition rates depend 
exponentially on distance and energy. The typical hop- 
ping distance and energy are functions of the temperature 
and the mechanism is known as variable-range hopping. 
The non-interacting version of this phenomenon was ex- 
plained forty years ago by Mott pQ. When interactions 
between electrons are taking into account the problem be- 
comes much more complicated. In the case of Coulomb 
interactions, Efros and Shklovskii (ES) proposed an ex- 
tension of the non-interacting case in which the constant 
density of states (DOS) of the non-interacting case is re- 
placed by the DOS of the interacting one [2J. This DOS 
presents a gap, known as the Coulomb gap, around the 
Fermi level [3j which limits the hopping of electrons at 
low energies. The shape of the Coulomb gap depends on 
the dimension of the system d and is given by: 



id— 1 



(1) 



for d > 1. In this formula the energy e is measured 
with respect to the Fermi energy. For d = 1 there exist 
logarithmic corrections to the DOS (1). One implicit ap- 
proximation performed in the ES argument is that weak 
quantum effects do not alter the shape of the DOS [5]. 

In this work we study the effects of a small transfer 
energy contribution t in the shape of the Coulomb gap. 
To generalize the definition of the classical Coulomb gap 
to quantum systems we employ the local DOS. The par- 
ticle contribution (e), and the hole contribution p^(e) 
to this local DOS at site i are [B]: 



P, 



(2) 



^(e)=^26(e + E^-Eo)\{^\ci\ih)\ , (3) 


where Eq, \i/)q) are the eigenenergy and eigenstate of the 
ground state in the grand canonical ensemble which we 
assume contains n electrons. The energy of an eigenstate 
a with an electron added to the ground states is denoted 



by and its corresponding wave function by \ipa)- F° r 
the energy and the wave function of an eigenstate (3 with 
n — 1 electrons we use Eg and \ipg) respectively. We 
notice that p\{e) and p^{e) have been defined in such a 
way that a positive energy is needed to add a particle 
while a negative one is necessary to subtract a particle. 
Using these particle and hole contributions, we express 
the average of the local DOS: 



P(e)-(p^e) + p1^)) 



(4) 



where the brackets mean average over different sites i or 
disorders. This expression reduces to the classical DOS 
when quantum effects are ignored. In the following, we 
refer to p(e) just as the DOS. 

An intuitive argument to describe the changes in the 
Coulomb gap due to quantum effects is given in Ref. [7|. 
In the classical case, the restriction on the DOS near the 
Fermi level is due to long-range interactions. If quantum 
effects only produce a small increase of the localization 
length, the nature of those long-range interactions should 
remain the same at scales larger than the localization 
length and so at energies (in units of e 2 ) smaller than 
its inverse. Then, one can conclude that the shape of the 
DOS close enough to the Fermi energy should not change. 
This picture agrees with the experimental results on tun- 
neling conductance |8l9j , although the relation between 
the DOS and the tunneling conductance is not yet fully 
understood [10] . We want to emphasize that this argu- 
ment does not seem to depend on the dimensionality of 
the system and the same behavior may be expected for 
any dimension. 

The first attempt made to include quantum effects on 
the Coulomb gap was performed by Vignale and cowork- 
ers pj]. They compute the DOS for non-interacting elec- 
trons. But, they chose an on-site energy distribution 
given by the classical Coulomb gap instead of being ran- 
domly distributed. Using the results obtained for this 
calculation, they claimed that the Coulomb gap is un- 
stable under the inclusion of arbitrarily small quantum 
effects in two-dimensional (2D) systems. In three dimen- 
sions, they found that the gap becomes narrower and 
quantum effects produce a non-zero density at the Fermi 
level. After that, Vignale P~2J used a Hartree-Fock based 
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approximation to obtain that the gap narrows for 2D and 
three-dimensional (3D) systems. In Ref. [13] a coherent 
potential approximation was employed to conclude that 
the gap narrows but the DOS at the Fermi energy re- 
mains negligibly small when weak quantum effects are 
considered in two and three dimensions. Schrciber et 
al. [7] studied the 3D case within a Hartree-Fock based 
method. Their results show that weak quantum effects 
narrow the gap and are compatible with a zero DOS at 
the Fermi energy. The narrowing of the gap with a finite 
t seems to be the natural tendency at least in 3D sys- 
tems, since the gap must be logarithmically small in the 
extended phase [Hj. However, Jeon et al. [TS] applied 
the same method to one-dimensional (ID) and 2D sys- 
tems and their results show that the slope of the classical 
gap near the Fermi level decreases upon including weak 
quantum effects. This is in contradiction with previous 
results |11)12|13] in two dimensions. We note the lack 
of agreement on whether weak quantum effects widen or 
narrow the Coulomb gap and on whether they produce 
or not a finite DOS at the Fermi level. Neither is clear 
if quantum corrections depend on the dimensionality of 
the system. 

We study the DOS given by Eq. Q for systems com- 
posed by many electrons when disorder and long-range 
interactions are important but quantum effects are weak. 
Our main goal is to study 2D systems with Coulomb 
interactions, but we also analyze ID systems with a 
v(r) = 1/y/r interaction. This interaction produces a 
linear gap in the classical DOS of ID systems, the same 
as the Coulomb interaction in the 2D case. This linear 
gap is easier to analyze numerically than the logarithmic 
gap produced by a Coulomb interaction in one dimen- 
sion. Furthermore, we will see that our results for ID 
systems show the same tendency that the ones found in 
Ref. [H] for Coulomb interactions. 

The structure of this paper is the following. In section 
|ri| we introduce the model and the basic numerical tools 
which shall be used. After that, in section [Til] we present 
results obtained with exact diagonalization. In section 
IV a perturbative analysis of the DOS is performed and 
the results of this analysis are used to compute the DOS 
for much larger system sizes than those obtained by exact 
diagonalization. Finally, our results are briefly summa- 
rized and discussed in section Ivl 



II. MODEL AND METHOD 



We study a spin-less fermionic system in a regular lat- 
tice with long-range interactions in which translational 
symmetry is broken by an on-site random potential and 
particles can tunnel to nearest neighbors sites. We use a 



tight-binding Hamiltonian: 
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where <pi is a random site energy with an uniform dis- 
tribution in the interval [— W/2, W/2], t is the hopping 
parameter, V is the strength of interaction, /i is the chem- 
ical potential and K is the compensation constant. The 
parameter a is chosen depending on the dimension as ex- 
plained in the introduction. That is, a = 1/2 is used in 
ID and a = 1 for 2D systems. In our computations we 
implement periodic boundary conditions in all the direc- 
tions and the interaction between two sites is set by the 
minimum-image convention. We take an interaction en- 
ergy V = 1, which fixes the energy scale, and we consider 
a disorder energy W = 2. We are mainly interested in 
a regime where the localization length is small compared 
with the system size, but non zero. We use a compensa- 
tion constant K and chemical potential /i depending on 
the filling to ensure the neutrality of the system. Most 
of our work has been performed at half filling for which 
K = 1/2 and /i — 0. For this value of K the DOS is sym- 
metric with respect to the Fermi energy. We shall take 
advantage of this property to only show half of the DOS 
when presenting results. We also compute in section |IV| 
the DOS at one third filling. In this case, we set K = 1/3 
and \x = —1/3. 

We analyze the two alternative approaches employed 
to obtain the DOS, Eqs. ^ and ((3). The first one is 
based on exact diagonalization. The largest system sizes 
considered with this method are L = 24 for one dimen- 
sion and N = 4 x 6 for two. In order to compute the 
DOS via exact diagonalization, we proceed as follows. 
We search for the ground state in the grand canonical 
ensemble which we assume contains n particles. As we 
will discuss later, this ensemble minimizes finite size ef- 
fects. We are interested in an energy region close enough 
to the Fermi level, so we only calculate the S smallest en- 
ergy levels for the cases corresponding to n + 1 and n — 1 
particles, with S ~ 50. Such a partial diagonalization is 
carried out using ARPACK open library [16] which are 
based on a variant of the Lanczos algorithm and allows to 
efficiently obtain a few eigenvectors and its eigenenergies. 
In Ref. [T7], we provide a thorough study of the validity 
of the method depending on the number of eigenstates 
considered. We run this procedure for at least 3000 real- 
izations of disorder to obtain an average of Eqs. ^ and 
([3]). Finally, by working in the grand canonical ensem- 
ble and choosing the minimum energy eigenstate for all 
the occupations, we ensure that the electron energies are 
always positive (hole energies negative) and this avoids 
an artificial filling of the DOS at the Fermi energy. An 
analogous problem in the classical limit is discussed in 
Ref. [IB]. 

In the second approach we use a numerical algorithm 
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FIG. 1. (Color online). The DOS as a function of energy for 
ID (left) and 2D (right) systems. The hopping is t — 0.2 and 
the sizes are specified in each panel. 



FIG. 2. (Color online). The DOS as a function of energy for 
size L = 20 in ID (left) and N = 4 x 4 in 2D (right) systems 
for several values of hopping as indicated in each panel. The 
thin black lines represent much larger sizes for the classical 
case and are shown as references. 



to obtain the classical ground state, i.e. the lowest energy 
configuration of the t = limit and a fairly complete set 
of low-energy configurations, needed to recalculate the 
new ground state in the presence of quantum effects as 
will be explained in section [TV] To compute these classi- 
cal states, we repeatedly start from states chosen at ran- 
dom and relax each sample using a local search procedure 
[19] . In an iterative process, we look for configurations of 
lower energies differing by one-electron or compact two- 
electron jumps and always accept the first of such states 
found. The procedure stops when no lower energy neigh- 
boring states exist, which insures stability with respect to 
all one-electron jumps and compact two-electron jumps. 
We then consider the set of metastable states found by 
the process just described and look for the sites which 
present the same occupation in all of them. These sites 
are assumed to be frozen, i.e. they are not allowed to 
change occupation, and the relaxation algorithm is now 
applied to the unfrozen sites. The whole procedure is re- 
peated until no new frozen sites are found with the set of 
metastable states considered. The lowest energy config- 
uration of the set of metastable states is consider to be 
the classical ground state. The set of low-energy states is 
completed by generating all the configurations that dif- 
fer by one or two electrons transitions from any of the 
metastable configuration state found before. 



III. EXACT DIAGONALIZATION 

As the sizes involved in the exact diagonalization are 
small, let us start with a finite size analysis to check 
that valid conclusions still can be drawn. In Fig. [T] we 
show the DOS for hopping parameter t = 0.2 in one 
(left) and two dimensions (right). The sizes employed are 



L = 12, 16, 20 in one dimension and TV = 4x3, 4x4, 4x 
6 in two dimensions. The symmetry of the DOS with 
respect to the Fermi energy at half occupation is taken 
into account and only half of the DOS is represented. In 
the classical case, finite size effects produces a filling of 
the DOS at the Fermi energy depending on the linear 
size of the system L as p (e = 0) ~ 1/L for both ID and 
2D systems [18]. As in this classical case, the results 
shown in Fig. [T] indicates that finite size effects for small 
hopping are important for energies close to the Fermi 
level. But they are not so important at intermediate 
energies inside the Coulomb gap. Indeed, we can observe 
in both panels appearing in Fig. [T] that there is a region 
around the middle of each branch of the gap where the 
changes in the DOS are small for the sizes considered. We 
take advantage of this fact to extract information about 
the behavior of an intermediate region of the gap when 
hopping is small but non-zero. 

We study now the dependence of the DOS on the hop- 
ping parameter. In the left panel of Fig. [2] the DOS 
as a function of energy is shown for several values of 
the hopping parameter for a ID system of size L = 20. 
From top to bottom, the different lines correspond to 
t = 0, 0.1, 0.2, 0.3 respectively. The thin line is the 
classical DOS for a size L = 1024 and it is included only 
as a reference. We can see that the gap in the DOS 
corresponding to an intermediate region widens for finite 
values of the hopping parameter. This is a surprising be- 
havior according to the intuitive picture explained in the 
introduction, but it is similar to the results found in ref- 
erence [15] for ID systems. Another remarkable feature 
is that the shape of the gap remains pretty much linear 
for the values of t shown. Near the Fermi level the gap is 
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not linear, but one can ascribe this to finite size effects. 

The right panel of Fig. [2] shows the DOS as a function 
of energy for t = 0, 0.1, 0.2, 0.3 in the 2D case for 
size N = 4 x 4. The thin black line corresponds to the 
classical DOS for a size N = 33 x 33 and it is included as a 
reference. Similarly to the ID case, we can observe that 
the gap in the DOS widens at an intermediate energy 
region as t increases up to t = 0.2, though the changes 
are smaller than in the ID problem. We stress that the 
tendency for t = 0.3 seems to be the opposite that the one 
for smaller hopping. That is, the DOS at t = 0.3 seems 
to narrow respect the DOS at t = 0.2. These results are 
again compatible with reference [15] but in contradiction 
with all the other previous works treating the 2D case 

mm- 



IV. PERTURBATION THEORY 

Our direct diagonalization method was constrained to 
small systems. In this section, we develop a novel ap- 
proach based on perturbation of the hopping parameter 
which will allow us to calculate the DOS for larger sizes 
and obtain reliable data for the DOS at Fermi level. Usu- 
ally, perturbation theory is difficult to apply to the kind 
of interacting system we are treating. In fact, a small per- 
turbation of a hopping term between two given sites af- 
fects particles which are arbitrarily far away of those sites 
due to the long-range character of the interaction. Then, 
the perturbative analysis must be performed over the full 
many-body space instead of the one-particle states. We 
have overcome these difficulties for small enough t by 
only considering the most relevant hopping terms when 
computing the local DOS at a given site i. 

In the classical case, t = 0, all the eigenstates have a 
well defined occupation in each site. As a consequence, 
there are well defined energies associated with each site 
of the lattice. To be more precise, the energy of site i is 
defined as: 




(6) 



where rii is the occupation of site i, Eq the energy of the 
ground state and Ef, E^ are the energies of the eigen- 
states with a particle or hole respectively created at site 
i to the ground state. The DOS for a classical system 
is just the density of these site energies. When quantum 
effects are introduced, this picture can be significantly 
altered. Indeed, the eigenstate of the Hamiltonian for 
non-zero hopping do not have a well defined occupation 
in each site and, in general, it is not possible to define a 
site energy. Nevertheless, we will see that a picture sim- 
ilar to the classical case still holds when the localization 
length is small enough. 

We compute the DOS treating the hopping term as 
a perturbation to the classical Hamiltonian. To do so, 
we can express the particle and hole DOS, Eqs. ^ and 



([3]), at site i as a power series of the hopping t. The 
first non-zero correction for these equations is of second 
order. The particle DOS at site i up to second order on 
the hopping parameter is: 



\P\ (*)](*) = t 2 J2He-K (0) + E (0)) ^ 

+ J26(e-E*(t) + E (<)) Ai(Q) + O (t 3 ) 

a 

(7) 



where A l a (t) 



are the matrix ele- 



(r a (t)\4\Mt)) 

ments appearing in Eq. ^ as a function of the hopping. 
The energy of the ground state up to second order cor- 
rection on t is E (t) and E? (t) is the energy of a state a 
with one more particle than the ground state up to sec- 
ond order corrections on t. An equivalent formula holds 
for the correction of the hole DOS Eq. ([3 1. 

At second order on t, an electron can hop at most to 
a nearest neighbor site. Roughly speaking, an electron 
hopping from site i + r to i + r ± I changes the classical 
energy of the site i as a dipolar contribution when r is 
large enough. Then, the changes in the classical energy 
at site i due to the hopping between sites i + r, i + 
r ± 1 should fall off faster than l/|r| 2 . We compute the 
DOS including only the dominant hopping terms, that 
is, the hopping terms between the site where the particle 
is created (destroyed) and its nearest neighbors. As we 
will see, this is an excellent approximation when hopping 
is small. 

The first non-zero correction to the matrix elements 
A l a due to nearest neighbors hopping is: 



Ai(t) = Ai(0)+t 2 J2 Ct a j + O(t 3 ), 



(8) 



where the sum in j runs over all nearest neighbors of site 
i and the C^ J are real constants. Using the standard 
perturbation theory it is possible to prove that these 
constants fulfill C^ 3 — —C 3 ^ 1 . Due to this property, the 
first term of the right hand side of Eq. ^ cancels when 
the particle DOS is averaged over all the sites in a sam- 
ple. We can also simplify the second term. Indeed, the 
matrix element A l a (0) is non-zero only if the eigenstate 
a corresponds to a state with an electron added at site 
i on top of the ground state. Then, our approximation 
leads to a DOS given by: 



p{e) = {8{e-e i (i))>, 



(9) 



where et(t) is the classical energy of site i, Eq. plus its 
second order correction due to nearest neighbors hopping 
and the brackets mean the average over sites. As in the 
classical case, each site contributes to the DOS in Eq. ^ 
with a well defined site energy. 

We have to determine the new site energies €i(t). In 
principle, the disorder breaks any degeneracy between 
the classical eigenstate and the standard-second order 
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FIG. 3. The DOS as a function of the energy for t = 0.1 
computed with the perturbation method (black solid line) and 
with exact diagonalization (red dashed line) for L — 16 in ID 
systems (left) and N — 4 x 4 in 2D systems (right). 
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FIG. 4. The DOS as a function of the energy for t = (black 
solid line) and t = 0.1 (red dashed line) for ID systems of size 
L = 256 (left) and 2D systems of size N = 10 x 10 (right) 
both of them computed with the perturbation method. 



perturbation theory could be employed to compute these 
energies. However, the existence of classical eigenstates 
with very similar energies can produce numerical insta- 
bilities. To avoid this problem, for each pair of nearest 
neighbors sites i and j we use an exact diagonalization 
scheme equivalent to the perturbation of a degenerate 
level. This hopping term only produces corrections in 
the classical eigenstates when n, 7^ Uj. Depending on 
the occupation of the ground state in these two sites, the 
hopping produces a correction either on the ground state 
or in a state with a particle or hole added on site i. For 
each of these cases, the correction can be computed by 
one of the eigenvalues of a two level operator Hij. A 
compact way of expressing Hij depending on the occu- 
pation and classical site energies of i and j is: 



Hij = {In, - 1) 







-t 



(10) 



where the factor (2rij — 1) is introduced to match the sign 
of the energy correction for particles (positive) and holes 
(negative) and 




express the site energy including all these corrections as: 



Ci(t) = £i + 



where the summation in j runs again over nearest neigh- 
bors of site i 

We notice that the ground state for non-zero hopping 
t can be different from the perturbed classical ground 
state. Indeed, some other low-energy state is likely to 
have the smallest energy when a finite hopping is taken 
into account. To deal with this problem we compute 
the quantum corrections of the first S classical states 
and we chose as the new ground state the state with 
the minimum total energy. We found that S = 10000 is 
enough for system sizes up to L ~ 300 in ID and up to 
N ~ 10 x 10 in 2D systems for the values of t considered. 
Then, the site energies given by Eq. ( 12 ) are computed 



using the d% t j factors associated to the zero order of this 
new ground state. 



if n 
if n 

— Ej + £j + 1 if Hi > Uj 



j t 1 if rii < rij 



(11) 



Finally, the site energy of % is the eigenvalue of Hij which 
reduces to e; when t = 0. 

We sum all the corrections given by Hij for j being a 
nearest neighbor of site i. This procedure can be viewed 
as a partial re-summation of the full perturbation serie in 
order to avoid numerical instabilities, although it is only 
exact up to the lowest non-zero order. Then, we can 



A. Results 

We begin checking the validity of our approximation by 
comparing with the exact DOS. In the left panel of Fig. 
[3j the DOS computed with exact diagonalization and the 
result from perturbation theory are shown as a function 
of the energy for ID systems for hopping t = 0.1 and 
size L — 16. At the right panel, the perturbed and exact 
DOS for 2D systems are presented for hopping t = 0.1 
and size N = 4 x 4. As we still be at half filling where 
the DOS is symmetric respect the Fermi energies, only 
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half of the DOS is shown as we did in the previous sec- 
tion. The agreement between both methods is remark- 
ably good in ID and 2D systems. For ID systems, we 
have checked that the perturbation method gives sensible 
results even for t — 0.2, while for 2D systems the results 
are not so good. We have also measured the localiza- 
tion length in ID systems obtaining £(f = 0.1) ~ 0.8 and 
£(i = 0.2) 1.4 [T7]. We conclude that the perturbation 
method employed is quite accurate for small values of £, 
of the order of the lattice spacing. In higher dimensions 
£ increases faster than in ID systems, thus, the range of 
validity of the approach is smaller. 

The DOS as a function of the energy for t — and t = 
0.1 is shown using ID systems of sizes L — 256 (left panel 
of Fig. |4| and in 2D systems of sizes N = 10 x 10 (right 
panel of Fig. Hj). First, as the sizes are much larger than 
the one used for the exact diagonalization method, the 
finite size effects are also much smaller. So now, we can 
extract a more reliable information about the behavior 
of the DOS near the Fermi energy. In Fig. gj we can 
appreciate that the DOS near this Fermi level does not 
seem to be affected by a small hopping in ID neither 
in 2D systems. Indeed, we have represented the DOS 
at zero energy, p (e = 0) , as a function of the inverse of 
the linear size, 1/L, for 2D system of linear sizes L — 
4, 6, 8, 10 finding that the results are compatible with a 
p (e = 0) = at the thermodynamical limit. On the other 
hand, we can appreciate in Fig. [4] that the DOS widens 
for a intermediate energy region as it was also found using 
exact diagonalization. Thus, the results indicate that the 
inclusion of a weak quantum perturbation reduces the 
slope of the Coulomb gap in ID and 2D systems. 

Finally, we check that our findings do not depend on 
the degree of filling. We compute the DOS at one third 
filling where the compensation constant is chosen as K — 
1/3 and chemical potential p = —1/3. In Fig. [5j we can 
see the DOS for hopping values t = (black solid line) 
and t = 0.1 (blue semi-dashed line) as a function of the 
energy calculated for a ID system of size L = 90 at one 
third filling and setting the Fermi level at the origin of 
energies. We also compare the DOS at this filling with 
the one that has already been calculated for half filling. 
The red dashed line is the DOS for t = 0.1 as a function 
of the energy at half filling for a system of size L = 90. 
We see that the DOS for hopping t = 0.1 at both fillings 
have the same shape near each of its respective Fermi 
energies. We have checked that the same result holds in 
2D. Then, the shape of the DOS near the Fermi energy 
for a small but non-zero hopping is independent of the 
compensation as it also occurs with the classical Coulomb 
gap [S]. 



V. CONCLUSIONS 

Our numerical results indicate that weak quantum ef- 
fects produce a widening of the Coulomb gap in the DOS. 
This widening is due to level repulsion. In the ground 
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FIG. 5. (Color online). The red dashed line represents the 
DOS as a function of the energy for t = 0.1 at half filling 
(p. — 0) for ID systems of size L — 90 computed with the 
perturbation method. The blue semi-dashed line is the DOS 
for t = 0.1 obtained also with the perturbation method and 
the black solid line is the classical DOS, t = 0, both of them 
computed for size L = 90 at one third filling where the com- 
pensation constant is if = 1/3 and the chemical potential 
/i = -1/3. 



state, when two nearest-neighbor sites have opposite oc- 
cupancy, the hopping term mixes two configurations: [om] 
and [»o]. Due to level repulsion the difference in energy 
must be larger than t. The energy of the ground state is 
always reduced by quantum corrections, while configura- 
tions with one more or one less particle are not affected 
(up to 0{t 2 ) ). This effect contributes to open the gap, 
as occupied sites will reduce its energy while empty ones 
will increase it. If two nearest sites are both empty in the 
ground state [oo], quantum corrections will affect only 
configurations with one more particle. Indeed, the site 
with lowest site energy will reduce it (due to level repul- 
sion) and the other one will increase it. So, for an empty 
site near the center of the gap this second effect will pro- 
duce, on average, a reduction of the site energy (the gap 
closes) as the chances are that the empty neighbor will 
have larger site energy. For ID and 2D we expect this 
second effect to be much smaller than the first one be- 
cause the probability of finding two empty sites with site 
energies within the gap and difference smaller than t is 
very small. The first effect does not require that both site 
energies should be inside the gap, so it should be more 
likely. In 3D systems, as the gap is parabolic, this second 
mechanism is more effective, so we expect the widening 
of the gap to be much smaller. 

The previous arguments are not rigorous as quantum 
effects, for a finite t, could change the configuration of 
the ground state. More properly we should use a self- 
consistent approximation taking into account the quan- 
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turn corrections. The net effect of this quantum correc- 
tions is the hardening of the constraints for the ground 
state. In particular an occupied site with energy e,-(t) 
and a nearest neighbor empty site with energy ej (t) must 
satisfy the condition 

ej (t)-ei(t)-l>t. (13) 

We may notice that the widening of the Coulomb gap 
with transfer energy t, specially on low dimensional sys- 
tems, should be taken into account in the interpretation 
of tunneling experiments in the strongly localized regime 
[819] . It may also affect variable range hopping conduc- 
tance. A calculation based on the classical DOS |20|21j 
will overestimate the conductance, since the widening of 
the gap will not be taken into account. 

In summary, we have studied quantum effects in the 
DOS of strongly localized regime with a l/\/r interac- 
tion in one dimension and a Coulomb interaction in two 
dimensions. Two different numerical methods have been 
employed. First, an exact diagonalization algorithm has 
allowed us to compute the DOS for small sizes. We have 
obtained that the DOS in one and two dimensions re- 
mains approximately linear for small but non-zero hop- 



ping. For ID systems, a monotonous reduction in the 
slope of the DOS near the Fermi level has been found 
when hopping increases from t=0to£ = 0.3. For 2D 
systems, the Coulomb gap also widens when the hopping 
is increased from t = to t = 0.2, but this tendency is 
reversed from there on. Using a novel perturbation ap- 
proach, we have been able to compute the DOS for much 
larger system sizes than for exact diagonalization using 
hopping t — 0.1. The results confirm that the slope of 
the Coulomb gap decreases when including weak quan- 
tum effects and that there is not an appreciable filling of 
the Coulomb gap at the Fermi energy. Finally, we have 
computed the DOS at 1/3 filling for ID systems. Our 
result indicates that the behavior of the DOS near the 
Fermi energy is independent of the degree of filling. 
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